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Abstract 


A nonperturbative analytic solution of the high charge and energy 
(HZE) Green s function is used to implement a computer code for 
laboratory ion beam transport in multilayered materials. The code 
is established to operate on the Langley nuclear fragmentation model 
used in engineering applications. Computational procedures are es- 
tablished to generate linear energy transfer (LET) distributions for a 
specified ion beam and target for comparison with experimental mea- 
surements. The code was found to be highly efficient and compared 
well with the perturbation approximation. 


Introduction 


Green sanction was identified as the likely means of generating efficient high charge and 
energy (HZE) shielding codes for space engineering that are capable of being validated in 
laboratory experiments (ref. 1). A derivation of Green’s function as a perturbation series 
gave promise for the development of a laboratory-validated engineering code (ref. 2) but 
computational inefficiency provided a major obstacle to the code development (ref 3) More 
recently, nonperturbative approximations to the HZE Green’s function have shown promise in 
providing an efficient validated engineering code (refs. 4 and 5). Described in the present report 

° rat0ry , C0de USmg a non P erturbative Green’s function to derive linear energy transfer 
(LET) spectra behind multilayered targets for ion beams with Z (charges) < 28 corresponding 
to the major components of the galactic cosmic ray spectrum. 


Green’s Function for a Single Medium 

We restrict our attention to the multicharged ions (atomic number Z 7 and atomic mass A<) 
tor which the Boltzmann equation may be reduced (ref. 6) to 3 

d d 


dx 


QE S j ( E ) + °j 


= J2 a jk 4>k(x,E) 

k 


( 1 ) 


where fy(x,E) is the ion flux at x with energy E (in MeV/amu), Sj(E) is the change in E 
per unit distance, Uj is the total macroscopic reaction cross section, and cr lk is the macroscopic 
cross section for the collision of an ion of type k to produce an ion of type j. The solution to 
equation (1) is to be found subject to the boundary condition 


4>j (0, E) = fj (E) ( 2 ) 

which, for laboratory beams, has only one value of j for which fj(E) is not zero and where ME) 
is described by a mean energy E 0 and energy spread a such that 


fj(E) = 


~7f=~ ex P 
V Zn a 


ziE_z e q ) 2 

2a 2 


( 3 ) 


The usual method of solution is to proceed toward solving equation (1) as a perturbation series 

(rets. 1 and 6). In practice, the computational requirements limit the usefulness of the technique 
for deep penetration (ref. 3). 



The Green’s function {G jm ) is introduced as a solution of 




Gjm (x, E, Eo) — ^ Vjlc Gkm ( x > E, Eo) 

k 


(4) 


(5) 


subject to the boundary condition 

G jm (0 ,E,E 0 ) = S jm 6 {E-E 0 ) 

where 6 jm is Kronecker’s 6 and 6 (E - E 0 ) is Dirac’s 6 function. The solution to equation (1) is 
given by superposition as 

( 6 ) 


<t,^x,E) = ^,Ja jk (x,E,E')f k (E')dE' 


If GUx E E') is known as an algebraic quantity, the evaluation of equation (6) may be 
accomplished by simple integration techniques, and then the associated errors m numencally 
solving equation (1) are avoided (ref. 7). 

The above equations can be simplified by transforming the energy into the residual range 

( r j) 38 E 

(7) 


= JdE'/SjiE') 


and defining new field functions as 

t l>j (x, Tj) = Sj {E) 4>j (x, E) 

Gjm (x, Tj, r'm) = $j ( E ) G J™ (x, E, E') 
fj (o) = Sj (E) fj (E) 

Thus, equation (4) becomes 

9 9 


+ Gjm {x, fj , r'm) = Gkm (*. r fc- r m) 

\dx dr 3 V k k 


with the boundary condition 

Gjm (0 ,Xj, r m) = tijmb ( r j — r m) 
and with the solution to the ion fields given by 


ip. 


oo 

(x, rj) = ^2 f Gjm (x, Tj , r'm) f m ( r m) rfr m 

in *L 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Note that Vj, which is the range scale factor as VjTj - i-Wm, is taken as Uj Zj/Aj. The 
solution to equation (11) is written as a perturbation series 

(14) 


Gjm (Xj r j i r m) ^ ' Gjm i x,r j ’ r m) 


2 


in which 


,( 0 ) 


and 


(*» r i> r m) = 9 U) 6jm t(x + rj - r' m ) 

(15) 

nG) ( t \ ~ v j G jm g(j,m) 

(16) 


where £?j*j(x, rj,r' m ) is zero unless 


( r j + x ) < r'm < —rj + x 


(17) 


for u m > uy If Vj > v vu as can happen in neutron removal, the negative of equation (16) is used 
and the upper and lower limits of equation (17) are switched. The higher terms are approximated 


as 


^fm ( x - r j> r m) » ^2 


‘ i a k i _ x m 9 U, *1 , *2> • ‘ • t *t-l, m ) 


x (% - Vj) 

g(j) = e-*i x 

9(jl,j2,---,jn,jn+l) = 1 1 1 ^n-\,jn) ~ 9 (J M2 . • • • » Jn-1, jn+l) 

^ n .+ 1 


In the above equations, 


and 


(18) 

(19) 

(20) 


(j) 

Note that Gj m (x,rj,r' m ) is purely dependent on x for i > 0, which we represent as t?-^(x). (See 
ref. 3.) In terms of the above, the solution to equation (1) becomes (from ref. 3) 


t j( x > r j ) = e ° jX h ( r j + x ) 


+ J2 Gjm ( x ) [^rn (r' ml ) - Fm (r' rnu )] 

m,i 


( 21 ) 


In equation (21), r mu and are given by the upper and lower limits of the inequality of 
equation (17). The symbol F m {r' m ) refers to the integral spectrum 


OO 



Fm ( r m) = 1 fm (r) dr 

(22) 


r' 

f rn 


We note that 

Fm { r in) — F m ( E ') 

(23) 

with 

oo 



F m (E') = J f m (E) dE 

(24) 


E f 
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and 


E f 


( 25 ) 


r' m = J dE/S m ( E ) 

0 

We now introduce nonperturbative terms for the summation in equation (21). 

First, we recall that the g function of n arguments was generated by the perturbation solution 
of the transport equation neglecting ionization energy loss (ref. 1) given by 


d 


9jm 


(#) — ^ ^ Ojfc 9km ( x ) 


subject to the boundary condition 
for which the solution is 


It is also true that 


9jm (0) “ fijrn 

9jrn i x ) = bjm 9 ( rn ) + a jm 9 O', m) + ■ ■ ■ 
9jm ( x ) = 9 i k ^ X ~ 9km ^ 


(26) 

(27) 

(28) 
(29) 


for any positive values of x and y. Equation (29) may be used to propagate the function g jm {x) 
over the solution space after which 


f \pi m ( a ') e J 

Qjm ( x , r J’ r rn) ~ e J ^ ( x + r ] r m ) + x (;/ m — jy^) 


(30) 


The approximate solution of equation (1) is then given by 
ipj (x, rj) w e~ a i x f {rj + x) 


+ E 


Vj 9jm ( x ) — e 


m ( r m«) ( r m 0j 


(31) 


x (r'm ^j) 
which is a relatively simple quantity (ref. 4). 

Green’s Function in a Shielded Medium 

The major simplification in the Green’s function method results from the fact that the scaled 
spectral distribution of secondary ions to a first approximation depends only on the depth of 
penetration as seen in equations (16), (18), and (30). Our first approach to a multilayered 
Green’s function will rely on this observation and assume its validity for multilayered shields. 

If we consider a domain labeled as “1” that is shielded by a second domain labeled as “2, 
the number of type j ions at depth x in 1 due to type m ions incident on domain 2 of thickness y 

9l2jm (z, y) = E 9l i k ^ 92km ^ ^ 

k 
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The leading term in equation (32) is the penetrating primaries as 

9Ujm (x,y) = e- 0 V x - a W6 jm + [g 12jm (x,y) - e~ a ' i x ~ a ^ v 6 jm 

in which all higher order terms are within the brackets of equation (33). 

The first term of the scaled Green’s function is then 

g i2jm ( x > r j > r 'm ) = e- a V x - a W '6 jm b [x + rj - (r' m - py)\ 

in which p is the range scale factor for the two media 


(33) 


(34) 


lj(E) 
R 2j ( E ) 


The ratio of range in water to range in aluminum for proton beams (eq. (35)) is shown 
in sketch A. We take a single value for p corresponding to 600 MeV/amu. The secondary 
contribution is similarly found by noting that equation (17) becomes 


( r j + ^ + PV ) < r rn — ~^~ r j + x + py 


(36) 


from which the average spectrum is evaluated. The full approximate Green’s function is then 


Q\2jm (x, V, rj , r' m ) « e °V X a ^ y 6 jm b{x + py + rj - r' m ) 

912 jm (x, y) - e~ a 'i x ~ a2 J y 6j r ,^ 


+ 


v 3 


{x + py) (v m - isj) 


(37) 


Equation (37) is our first approximation to the Green’s function in a shielded medium (two 
layers) and is easily modified to multilayers. We now consider the first spectral modification. 



E, MeV 
Sketch A 
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Showing that the first collision term has the properties 


^12jm rm ) ^ 


I "m - Vj\ 

j x ~ a 2jy 

y j a 2jm e 3 J 


{ r m T mu) 
( r m ~ r m £ ) J 




(38) 


r m - Vj I 

is easy; we use these properties to derive a simple correction for the average spectrum as 

„(D 


Q { S>. ( x , y, r 3 , r' m ) = - - gl2 J" t(X ’ y -~ + b jm (x, y) {r' m - f m ) 

* 12 jm\ >»> J' mj (x + py)( Vm-Vj) 

where ffi2j m ( x ’y) is the first collision term of e 9 uation ( 39 ) and 


(39) 


r 1 4- r' 

-l 1 mu ' ml 

r m~ 2 


(40) 


is the midpoint of f' m between its limits given by equation (36). The b jm term of equation (39) 
has the property that 

J t ^ b jm (x, y) (r' - r m ) dr' = 0 ( 41 ) 


thus ensuring that the first term of equation (39) is, indeed, the average spectrum as required. 
The spectral slope parameter is found to be 


UjUm (o Xjm e-^ x -^ y ~ o 2jm e 
b 3 m (X,2/) = (x + py) {vm - Uj) | r'm - Vj 


(42) 


A similarly simple spectral correction can be made to the higher order terms. The spectral 
correction given in equation (42) will be included in the present Green’s function code. 

LET Spectra for Laboratory Beams 

We use the boundary condition appropriate for laboratory beams given by equation (3). The 
cumulative spectrum is given by 


Fj W = o 


1 — erf 


(E-E 0 \ 1 

V n/2<t ). 


(43) 


The cumulative energy moment needed to evaluate the spectral correction is 


JE-E 0 X 

a 

(- E-Eof 

L 1 " 61 ( v/2a )J 

+ r—- exp 

2a 2 


Ej (E) = \e o 


The average energy on any subinterval {E\, E 2 ) is then 

Ej (Ei) - Ej (E 2 ) 
Fj {Ei) - Fj (E 2 ) 


(44) 


( 45 ) 
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The beam-generated flux is 


*3 (*, V, rj) = e a V x jV f. ( rj +x + py ) 

+ £sj£ (*,!/) [ p ™(r' m u) -Fm(r' ml j\ 

m,i J 

bj m ( x i V ) [ r m (E) — r m ] \F m {r' mu ) — F m {r' m ^\ (46) 

m L J 


w here is evaluated using equation (45) with E 1 and E 2 as the lower and upper limits associated 
with r mf and r' mu , respectively. 


The differential fluence spectra for a 600-MeV/amu 56 Fe beam with a 2.5-MeV/amu standard 
deviation incident on a water slab are shown in figure 1. A single layer of 5 cm of water is shown 
in figure 1(a) and a double layer of 2.5 cm of water followed by 2.5 cm of water is shown in 
figure l b) A consistency check is performed on the multilayered code by comparing it with 
e single-layered computation when the two layers are of the same size and composition The 
ratio of multilayered results to single-layered results differs by less than ±1 percent. The LET 
distribution is found by using the methods of reference 8. The corresponding LET spectra are 
shown in figure 2. The highest LET peak is due to the primary beam and the iron fragments. 

e successive peaks below iron are due to lower atomic weight fragments. The lowest LET peak 
consists of relativistic charge fragments of p, d, and t particles that are produced in abundance 
m HZE collisions (ref. 9). The peak near 10 MeV/cm consists of he- and o-particles that are 
also produced abundantly. Many of the HZE fragments are produced with a charge near the 
projectile charge, as Shinn, Townsend, and Wilson found earlier for hydrogenic targets (ref. 10). 
Note that no distinguishable difference exists between the LET spectra of the single- layered code 
(fig. 2(a)) and the multilayered code (fig. 2(b)) at the same penetration depth in water, which 
further demonstrates code consistency. 


series of evaluations for a (2.24-g/cm 2 ) lead-scattering foil is shown in figure 3. The 
ead-scattermg foil is usually part of the accelerator beam line (at the Lawrence Berkeley 
Laboratory Bevalac accelerator) with the result that the fragments from the lead target are 
seen as contamination. Clearly, these fragments must be modeled to properly interpret the 
attenuation of the beam in the water target in actual experiments. The corresponding LET 
spectra at various depths of a water target are shown in figure 4. The importance of the 
ragmentation in the scattering foil is seen in comparing figure 4(a) with figure 4(b). Note 
that fewer of the fragmentations result in fragments near the beam charge for the lead foil in 
comparison with the water target, as seen by comparing the relative magnitudes of the three 
highest LET peaks in figures 2(a) and 4(a). These differences are part of the reason why 
ydrogemc targets are important to space radiation protection. At greater depths the LET 
distributions begin to overlap, and distinguishing the different charge groups becomes more 
difficult. Such LET spectra will be compared with experimental measurements in the near 


Concluding Remarks 

A formalism for the evaluation of Green’s function in multilayered target configurations 
has been derived and a computer code generated. The code satisfies a consistency check for 
multilayered-material calculations when the layers are of uniform composition by comparing the 
results with the single-layered code. The importance of the multilayered code in the analysis of 
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experimental ion beams has been shown by demonstrating the effects of a lead-scattering foil in 
the Bevalac beam line. An analysis of such experiments is in progress. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
September 10, 1093 
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(b) 2.5 cm of H 2 0 followed by 2.5 cm of H 2 0. 
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Flg r 3 ' ?!ff, ent / ial 5'‘ enCe f ° r 525_MeV / amu 56Fe be am with 2.5-MeV/amu standard deviation after passing 
through 2.24-g/cnr lead-scattering foil and water target. 
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LET, MeV/cm 

(b) 5 cm of H 2 0. 


[eV/amu 56 Fe beam with 2.5-MeV/amu standard deviation after passing 
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